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In this paper initial state radiation corrections to double production of charmonium mesons on 
one-photon electron-positron annihilation at center of mass energy yfs = 10.6 GeV are studied. It 
is shown, that these corrections have noticeable effect and must be taken into consideration. 
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I. INTRODUCTION 

Since the discovery of the so called charmonium mesons 
(that is mesons consisting of c- and c-quarks, for example 
J/ip, ip{2S) or rjc) in 1974 [lj, |2| they play a significant 
role in our understanding of quantum chromodynamics 
(QCD). Production of charmonium requires creation of 
heavy cc pair with the energy greater than 2m c where 
QCD coupling constant is small enough to use pertur- 
bation theory. However, the subsequent hadronization 
probes much smaller mass scales of order m c v , where 
v is typical velocity of c quark in the charmonium rest 
frame. For J/ip, ra c v 2 is numerically of order Aqcd, so 
the production processes are sensitive to nonperturba- 
tive physics as well. Charmonium mesons are also very 
interesting from the experimental point of view, since 
the narrow resonances can easily be separated from the 
background. 

In 1986 Caswell and Lepage have proposed the 
Non-Relativistic Quantum Chromodynamics model 
(NRQCD) [10. This model exploits the fact that the 
velocity v of the heavy quark in the charmonium rest 
frame is small in comparison with the speed of light and 
all physical quantities can be expressed as the series in 
this small velocity and electromagnetic and strong cou- 
pling constants a and a s . Matrix elements of the four- 
fermion operators, used in this theory, can be determined 
phenomenologically. Both inclusive and exclusive pro- 
duction of charmonium mesons and their decays were 
studied using this theory U H Q, H El EH 

New difficulties have arisen as a result of recent mea- 
surements of inclusive charmonium production in e + e~ 
collisions by Belle and B ABAR collaborations [L?! Il3| . 
Their results were about an order of magnitude higher, 
than the predictions based on NRQCD or other models 

S El El El 113 . 

Another difference has arisen in studying the exclusive 
production of pair of charmonium mesons. For example, 
the NRQCD results for J/ipr) c production is 
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J/iprjc) = 2.31 fb 



and Belle results E3 are an order of magnitude higher. 
One of the pos sible explanations of this difference was 
proposed in [19j . The authors assumed that some of the 
Belle's J ji>r\ c signals could actually be the double pro- 
duction of J/'i/j meson with subsequent decay J/ij) — > r/cj 
and presented the value 



2 7 * -> 2J/ip) = 8.7 fb. 



Later in [3, |2(j it was shown that due to relativistic and 
higher order QCD corrections the cross section should be 
about 4 times smaller and the question of the difference 
between theory and experiment remains open. 

In this paper I propose that this difference can be ex- 
plained by the initial state radiation corrections. Indeed, 
the naive estimate of the effect of these corrections gives 
us the result 



a = cree 
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and the suppression caused by additional factor a will 
be compensated by large logarithm log sjm e . Another 
reason is that cross section of the reaction e + e~ — * J/tprjc 
as a function of center of mass energy y/s has a narrow 
peak near yfs ~ 6.5 GeV, where a (e + e~ — > J /ipr) c ) ~ 
30 fb and by emitting a hard photon we can return to 
this region. 

The rest of this paper is organized as follows. In the 
next section I give the diagrams and matrix element 
for the process under consideration. In section IIIII the 
method of calculation of the total cross section is de- 
scribed. Finally, distributions over vector charmonium 
energy and scattering angle are presented. 



II. MATRIX ELEMENT 

The diagrams for the process e~ (fci) e + (^2) — > 
J/ip (pi, E) j] c (j>2) 7 (k, e) are shown on figure^ The ver- 
tex of the transition of virtual photon with momentum 
p and its square s = p 2 into pair of charmonium states, 
denoted here by filled circle, was calculated in || ■ Some 
of the diagrams of this transition are shown on figure |3 
the others can be obtained from them by permutations. 
One can also add electromagnetic diagrams by replac- 
ing the gluon line in figure |2I a) by photon, but they are 
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suppressed by the factor a/a s and will not be considered 
here. On the other hand, as it was shown in 8] the purely 
electromagnetic resonance diagrams (for example, the di- 
agram^b)) must be included. The reason for this is that 
the additional factor a is compensated by r 2 = 4M 2 /s, 
where M = 2m c is the mass of the charmonium meson 
(the difference between J/ip and r\ c masses is neglected) 
and m c is c-quark mass. 

The analytical expression for the "fJ/ipr] c vertex was 
found to be 

J n = (J/^(pi,£)Vc{p2) |c7m c I°> = 

where e is the polarization vector of the vector charmo- 
nium J/tp, the coefficient A is 



A 
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N c = 3 is the number of colors, e c = 2/3 is the charge 
of the c quark (in the units of the elementary charge e) , 
and the constants 

= (Oi)^, o v = (o 1 ) Vc 

are the matrix elements of the NRQCD probability fac- 
tors. These constants can be determined from the exper- 
imental data on J ftp and r\ c decay widths @: 
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Using the table data on particle decay widths |2JJ we get 
the numerical values of probability factors, that will be 
used in this paper: 



0.335 GeV 3 



O n = 0.297 GeV 3 . 



The matrix element corresponding to the diagrams 
shown on figure ^ is 



M 



1 



(1) 



where k\ and k 2 are electron and positron momenta (the 
mass of the electron m e is neglected wherever it is possi- 
ble), k is the momentum of the final state photon, e M is 
its polarization vector, q a = k\ — k and qi, = k — k 2 are 
the momentum transfer in this two diagrams. The cross 
section found from and equals to 
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TT^fcifca) / M3 ( fc i + fc 2;Pi,P2,A:)^|X| 2 , 



where summation is performed over polarizations of ini- 
tial and final particles and for the Lorentz- invariant phase 
space the notation 



n \ n 



d$ n (P; Pl ,...,p n ) = 5[P-J2p* ( 2 ) 

\ i=l / i=l 1 

is used. 

III. TOTAL CROSS SECTION 

The Lorentz-invariant phase space J2J satisfies the re- 
cursion relation 

d<fr 3 (ki + k 2 ;pi,p2,k) = d'sd^ 2 {k\ + k 2 ;p,k) 
x d$ 2 (s;pi,p 2 ) , 

where p = Pi +p 2 is the momentum of the virtual photon, 
s = p 2 = s — 2y/sEk, and Ek is the final state photon 
energy in the laboratory frame. Leaving only the inte- 
gration over we get the differential cross section 



da _ 327687r 2 a 3 O^O,, 4M^ 



ds 19683 



(2as-4M 2 (a + 2a s )) 2 
x j — x 

s 4 

^s 2 + s 2 arctanh B e 
x I 2 — — — — - s + s 



From this equation one can easily notice two important 
items: 

• the cross section is enhanced by the factor 
arctanh j3 e ~ log s/m 2 as it was states in the in- 
troduction, 

• total cross section actually diverges in the limit 
s « s. This fact is the familiar infrared divergence 
caused by the emission of the soft photon, where 
one can not use the perturbation theory. To avoid 
this divergence one must put the cutoff S t on the 
energy of the photon: 

E k > 5 e . 

The plot for differential cross section da/dEk is pre- 
sented on figure 01 As stated above, in the limit of small 
photon energy it grows to infinity and one should put the 
cutoff. The bump in the region Ek~?> GeV corresponds 
to y/§ sj 6.5 GeV, where the cross section of the process 
7* — > J/ipr] c is maximal. 

Values of total cross section for some cutoff photon 
energies are: 

a w 2.84 fb for 8 e = 1 MeV, 
a « 2.234 fb for 5 e = 10 MeV, 
a « 1.72 fb for 8. = 100 MeV. 
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FIG. 1: Diagrams for e + e — > J/i/jr/c-y. 
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FIG. 2: Some diagrams for the transition of a virtual photon into charmonium pair. Other diagrams can be obtained from 
theses by permutations. 
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FIG. 3: Distribution over final state photon energy 

It is easily seen, that by the order of magnitude these 
cross sections equals to the cross section of the base pro- 
cess e + e~ — * 7* — > J/ijjric, so the initial state radiation 
corrections should be taken into account. 



Vk,X,2 



2M 

COS 6 k ,1, 2, 



1 _ T 2 

x x fc,l,2 ' 



where E\ , E% and E k are the energies of vector and pseu- 
doscalar charmonia and photon in the center of mass 
frame respectively; 9\, Qi and 9 k are the angles between 
the momentum of the initial electron and the momenta 
of these final particles. These variables take the values 
in the regions 

xi,x k e [-l;l], 

z k € [z™; z k ], 
where minimum and maximum values of z k are 



IV. DIFFERENTIAL CROSS SECTION 



z m,M _ 2(1 — Zi) 

k 2- zi\l- (3{xix k ±yiy k )\ 



For the calculation of the differential cross section over 

the vector charmonium energy it is convenient to intro- . . 

j . , .I, j. , c i j In the terms of these variables the squared matrix de- 

duce dimensionless variables according to tormula s T ... , , „ 

ment and Lorentz-mvariant phase space take the form 

z k,l,2 — 7= — 3 
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]M f . 134^0^ (3a,.- 0(3-3* - 1 4 _ , _ ^ 
6561/i 4 s 4 z£(l - z fc ) 5 (l - /3 2 .x£) [ 

4 (4 — /1 2 — 3 (3 + /9xia: fc ) 21 + (l + P 2 x\) z\) z k + 

(2 (5 - 2 Zl (1 - 0x lXk ) + 4) - (1 + 4)) - 




where /3 = \/l — M 2 / E 2 — y/\ — H 2 1 z\ is the velocity 
of the vector charmonium meson and /3 e = yl — 4to 2 /s 
is the velocity of the initial electron in the laboratory 
frame (this is the only place where I have to leave the 
non- vanishing electron mass). 

For integration of the differential cross section over z k 
it is possible to use the residue technique. The only 
d 4 a J ' dx\dz\dx z dz k poles are z k = 0, 1 in , so we can write 



d 3 a 



dx\dz\dxk 



E 

a* =0,1 



do 



res 



dxidzidxkdzk 



The residue at Zk = 1 is too lengthy to be presented here, 
but the residue at z k = 0, that plays most important 
role for z\ ~ 1 (that is in the limit of maximal vector 
charmonium energy) and X\ ±1 it is rather compact. 
In this limit we have 



B(zi) 



2zi + z{ 



The integration over all other variables was performed 
numerically. 

The distribution over vector charmonium energy is 
shown on figure 21 I n the limit of maximal vector char- 
monium energy (z\ rs 1) this curve grows to infinity cor- 
responding to analogous grows int the distribution over 
photon energy shown on figure |3J 

The distributions over scattering angle for different val- 
ues of S e are shown on figure This form of distribution 
is usual for color-singlet charmonium production. 



d 2 a 
dz\dx\ 



S( Zl ) [1 + B( Zl )x 2 ] 



n . , 2621447r 2 a 3 /3 CLO,, arctanh/3 e 
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(2 - n 2 - 2 Zl + z 2 ){3a sf i 2 - a(3 - ^ 2 )) 2 
l-zi 
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FIG. 4: Distribution over vector charmonium energy. 
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FIG. 5: Distribution over scattering angle for 5 e — 1 MeV 
(solid line), 5 e = 10 MeV (dashed line), and 5 t = 100 MeV 
(dotted line) 



